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PREFACE 


Geodetic satellites such as GEOSAT, SPOT, ERS-1, and TOPEX/Poseidon require 
accurate orbital computations to support the scientific data they collect. Until recently, 
gravity field mismodeling was the major source of error in precise orbit definition. How- 
ever, with recent improvements in these models, accurate modeling of the nonconser- 
vative forces has become a significant concern. In fact, the TOPEX/Poseidon mission 
requirements dictate that the mismodeling of the nonconservative forces of solar radia- 
tion, Earth albedo and infrared re-radiation, and spacecraft thermal imbalances produce 
in combination no more than a 6-cm radial rms error over a 10-day period. This re- 
quires the development of nonconservative force models which take the satellite's 
complex geometry, attitude, and surface properties into account. Specifically, for 
TOPEX/Poseidon, a "box-wing" satellite form has been investigated. This approach 
models the satellite as the combination of flat plates arranged in the shape of a box and 
a connected solar array. The nonconservative forces acting on each of the eight sur- 
faces are computed independently, yielding vector accelerations which are summed to 
compute the total aggregate effect on the satellite center-of-mass. Parameters associ- 
ated with each flat plate are adjusted to obtain a better representation of the satellite 
acceleration history. In order to test the validity of this concept, "micro-models" based 
on finite element analysis of the TOPEX/Poseidon spacecraft have been used to gener- 
ate acceleration histories in a wide variety of orbit orientations. These accelerations 
profiles are compared to those from the "box-wing" model. Tests indicate that modeling 
solar radiation pressure acceleration is relatively straightforward. However, the thermal 
imbalance modeling is made much more complicated given the satellite's complex 
attitude control law and its relation to the predicted temperature history for each surface. 
The results of these simulations and their implication on the ability to precisely model 
the TOPEX/Poseidon orbit will be discussed. 
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INTRODUCTION 


Mission/Science Overview 

The seasonal variations of the world’s ocean circulations, and how they influence the 
Earth’s climate will soon be investigated with a high-precision Earth orbiting altimetric 
satellite. This project, referred to as the Ocean TOPography EXperiment/Poseidon Mission 
(T/P), is a joint venture between the U.S.’s National Aeronautics and Space Administration 
(NASA) and the French Centre National d’Etudes Spatiales (CNES). The TOPEX/Poseidon 
spacecraft will be equipped with two radar altimeters (1 US, and 1 French), which will 
measure the ocean surface topography. 

By measuring the height of a satellite above the ocean surface with the use of a radar 
altimeter and subtracting the satellite’s height in geocentric coordinates, the sea surface 
topography in geocentric coordinates is calculated. Subtracting the geoid height and 
accounting for Earth/ocean tides yields a measure of the sub-satellite ocean dynamic 
topography. Knowledge of the dynamic topography over time is very important in that it 
consists mainly of surface geostrophic currents and ocean thermal responses which drive 
global weather patterns and their changes. 

The T/P spacecraft will orbit the Earth at an altitude of 1336 km, inclination of 66° and with 
nearly zero eccentricity. The period of the orbit is 1.87 hours and its groundtrace will 
repeat every 10 days in a "frozen" orbit. T/P is expected to be launched in July of 1992 
using the European Space Agency’s Ariane launch vehicle. 

Precision Orbit Determination 

Since the orbit of T/P will provide the reference frame from which the altimeter 
measurements are taken, any error in determining the satellite’s position will affect the 
accuracy of the direct interpretation of the altimetric measurements. In order to obtain 
direct measurements of dynamic topography to the degree of accuracy that is required 
by the scientific community for oceanographic studies, the orbit of T/P must be 
determined over contiguous 10 day periods to within 13 cm RMS accuracy in the radial 
component with less than 5 cm RMS geographically correlated error [Stewart et ai, 
1986]. Orbit determination of this accuracy has never before been undertaken or 
achieved for a satellite at T/P’s altitude. 

The common method for computing the radiation pressure upon orbiting satellites with 
the orbit determination software has been to ignore these rotating, attitude controlled, 
geometrically complex shapes and treat the satellite form as a symmetrically perfect and 
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rotationalty invariant sphere, or so-called "cannonball." A typical cannonball radiation- 
pressure model had a constant projected area with respect to the radiation source, and 
empirical scaling factors were commonly adjusted to improve the orbital fit to the tracking 
data. In addition, the thermal emission of the spacecraft itself was ignored. These 
approaches are not adequate to meet the required 6 cm radial rms error budget for 
modeling the radiation forces on T/P over a 1 0 day period. After considerable analysis 
of all surface force contributions, resultant models to be used in T/P orbit determination 
have been derived and are presented herein. 


Radiation Forces 

At the 1336 km altitude, solar radiation is the greatest nonconservative force acting on 
the spacecraft. In addition to direct solar radiation from the Sun, the Earth’s albedo and 
infrared (IR) emissions, along with infrared emissions from the spacecraft itself will perturb 
the orbital motion. These effects are described below. 

Solar Radiation 

The Sun emits a nearly constant amount of photons per unit of time, varying less than 
0.2% [Kivelson, 1986]. At a mean distance of 1 A.U., this time rate of flow of radiant 
energy per unit area, known as the solar constant or solar irradiance [Wilson and 
Hudson, 1988], is 


Gq = 1367 .7 W/m 2 


Since the Earth’s orbit is eccentric, the solar radiation flux received by the Earth varies 
by ±40.6 W/m 2 throughout its orbit around the Sun. The force produced by this radiation 
is by far the largest of the radiative effects. Also, for T/P, it is the largest nongravitational 
force acting on the satellite. 

Earth Albedo and Infrared Radiation 

A fraction of the solar radiation is instantly reflected off the land-water-snow-ice- 
atmosphere system as shortwave energy and the remainder is absorbed by the Earth. 
This absorbed energy is later emitted as longwave infrared radiation. 

The development given by Knocke et al., 1988 for modeling the radiation force due to 
Earth albedo and infrared emission has been adopted for this study. For the case of 
albedo, the amount of radiation received by a satellite due to an elemental Earth surface 
area dA&, will be, 
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where 6 0 is the angle between the unit normal of the elemental Earth surface area and the 
direction to the Sun, y is the angle between the Earth surface unit normal and the 
direction to the satellite, r is the distance from the Earth surface element to the satellite, 
and a is the albedo coefficient of the surface element. This model assumes that the Earth 
reflects and emits most of its energy in a diffuse manner and, therefore, each elemental 
Earth area radiates constant flux. In general, the a varies over the entire irradiated portion 
of the Earth. For this application, a longitudinally averaged model of albedo, as specified 
by Knocke et al. [1988] has been used. The average value of the albedo coefficient is 
0.34, which means the reflecting Earth surface elements radiate, on average, 465 W/m 2 . 

The infrared radiant energy flux density received by the spacecraft from an elemental 
Earth area dA#, at a distance, r, is 

G© ■ -^ cos YcM® 

* 47ir z 


where € is the emissivity of the surface element. As with the albedo, a simple latitudinal 
model is used for the emissivities. The average emissivity of the Earth is 0.68, which 
results in the surface elements radiating, on average, 232 W/m 2 . 

Spacecraft Radiation 

In addition to the radiation forces produced by the direct effects of incident radiation, a 
heated surface produces a thermal force due to its emission of infrared energy. A 
spacecraft surface with a temperature of T will emit longwave radiation at the level 

G*/c = eoTi 

where a is the Stefan-Bottzmann constant (5.670 x 10* W/m 2 K 4 ) and € is the emissivity 
of the surface. As the temperature of the surface changes throughout the orbit 
(especially due to eclipsing) the magnitude of the thermal force can change substantially. 
The surface itself changes because of the hostile orbital environment, and these effects 
must also be considered. 
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TOPEX/Poseidon Spacecraft 


The T/P spacecraft will fly in a near-circular orbit at an altitude of 1336 km and an 
inclination of 66 degrees. The eccentricity is controlled to maintain the argument of 
perigee at 90 degrees, producing a "frozen" orbital mission. This frozen orbit configuration 
minimizes height variations (with respect to the ellipsoid) and is very favorable for 
altimetric missions. Table 1 summarizes the orbital characteristics. 


Table 1. Orbital Parameters 


| Semi-major axis 

a 

7714 km 

| Eccentricity 

e 

0.00 

Inclination 

i 

66° 

Argument of Perigee 

6> 

90° 

Altitude 

h 

1336.00 km 

Nodal Precession Rate 

n 

-2.31 deg/day 

Period 

P 

1 1 2 min H 


The spacecraft consists of an instrument module, propulsion module and a single solar 
array. A representation of the spacecraft is shown in Figure la and the corresponding 
approximate projected areas of the spacecraft are summarized in Table 2. Note that there 
is a significant variation in the projected area-to-mass ratio depending on the view 
direction and that the spacecraft exhibits considerable asymmetries. 


Table 2. Projected Area Properties 


H 

X projected s/c body area 

A 

4.70 m 2 

Y projected s/c body area 

A 

8.18 m 2 

Z projected s/c body area 

A 

8.30 m 2 

| Solar array area 

A 

25.5 m 2 

| Spacecraft mass 

m 

2500.0kg 

g Max. proj. area-to-mass ratio 


0.140 cm 2 /g 

| Min. proj. area-to-mass ratio 


0.01 9 cm 2 /g 
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Attitude/Yaw Steering 

The TOPEX/Poseidon spacecraft has an intricate attitude control system because of its 
large, single-axis gimbal solar array. Perfect solar pointing (Sun incidence vector normal 
to the solar array) requires the spacecraft to yaw about its Earth-pointing Z-axis at rates 
that exceed the capabilities of the attitude control system. Therefore, Fairchild Space 
Company (FSC) designed a sinusoidal yaw command that achieves near-perfect solar 
pointing while remaining within attitude control system limits [Perrygo, 1987]. The 
algorithm is based on the coordinate systems shown in Figures 1 b and 2. The spacecraft 
body-fixed system origin is within the vehicle body, although not at the center of mass, 
with the positive Y-axis pointed opposite of the solar array axis, the positive Z-axis 
directed to the Earth nadir, and the positive X-axis orthogonal to the Y and Z axes to 
complete the right-handed system. The inertial system is centered at the geocenter with 
the X-axis normal to the satellite orbit plane, the Z-axis pointed in the direction of the 
Sun as projected into the orbit plane, and the Y 0 -axis normal to these axes. As shown 
in Figure 2, B’ refers to the angle between the Sun vector and the orbit plane and a is 
the orbit angle, measured from the Y c -axis. 

The actual yaw angle, ijr, of the spacecraft is rotated positively from the Xp-axis (the along- 
track direction) about the Z-axis and is determined from the B’ and Q. The multiple yaw 
steering algorithms the satellite follows during its orbit are summarized in Table 3 
[Zimbelman, 1989]. The B’ and n yaw steering mode boundary values specified are the 
nominal values selected for the mission and may be changed at any time by ground 
controllers. 
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Table 3. TOPEX/Poseidon Yaw Steering Modes 


Yaw Steering 
{ Mode 

S’ Region 

n Region 

Yaw Algorithm | 

Fixed 

+0. 1° SB's +15° 

All 

0 

O 

II 

&< 

-15° sB's-O . 1° 

All 

W = 180° 

B^ + 80° 

All 

II 

O 

o 


w 

(A 

i 

o 

o 

All 

o 

O 

<? i 
1 

II 

Sinusoidal 

B 7 > +15° 

All 

T=90° + (90°-B / )COSQ 

B'< -15° 

All 

Y = -9 0° - ( 9 0 0 +B') cos Q 

Ramp Up 

B'= + 15° 

90°iQ^180° 

Y=B'cos 2 Q 

B'=-15° 

27 0° £Qs360° 

T = -B / cos 2 Q-180° 

Ramp Down 

B /= +15° 

180° 0° 

Y=B'cos 2 Q 

B /= -15° 

1 

0°i£5i90 o 

T=-B / cos 2 Q-180° 

— — 


0°^B / <0.1° -180°^Q^0° Y = -90° [1+cosQ] 












































As Table 3 suggests, the yaw angle computation algorithm is extremely intricate at the 
yaw regime boundaries. When the on-board computer indicates that the spacecraft has 
crossed a B’ limit, it sets a flag to initiate the steering mode when it next crosses into the 
appropriate orbit angle transition regime. If the satellite crosses a B’ limit while in its 
corresponding orbit angle transition area, it will wait an entire revolution to commence 
yaw mode changes. 

The solar array pitch angle is defined as, 


0°<B'<15° 

-15 o <B'<0° 

P1js15° 


y=q 

y=h-q 
Y=arcta 


ri 


sinQ cosP / 


\cosQ cosP' cosY-sinY sin(J 
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where the pitch angle, y. rotates positively from the spacecraft X-axis about the Y-axis to 
orient the cells toward the Sun for optimum Sun pointing. 


Model Development 

The first step in a detailed analysis of the radiation forces acting on T/P was to accurately 
compute the radiation forces due to the Sun, Earth albedo, Earth infrared, and spacecraft 
thermal emissions upon T/P with the use of a finite element model of the spacecraft. This 
investigation generated what is referred to as the "micro-models"; a thorough explanation 
of this modeling effort is given by Rosborough and Antreasian [1991]. Since a precise 
thermal and radiative model of a spacecraft is necessarily computationally intensive, this 
micro-model, which served as a "truth" model, was computed offline. A relatively simple 
and less computationally intensive model, or macro-model, more suitable for precision 
orbit computations, was devised and tested to emulate the micro-model accelerations. 
This model was introduced in Marshall et ai, 1991. A graphical representation of this 
development is shown in Figure 1 . 

Micro-Models 

A radiation force is produced on a spacecraft surface if that spacecraft surface absorbs, 
reflects and/or emits radiant energy. Hence, forces will be generated by incident radiant 
energy from the Sun or Earth falling on a surface and by radiant energy emission from 
the surface. The detailed evaluation of the radiation forces is dependent upon accurately 
defining the flux that is incident on the various spacecraft surfaces and the resulting 
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interaction (absorption, reflection). The evaluation of these processes has been carried 
out using industry standard software by correctly orienting the spacecraft for a given 
orbital location and determining the incident fluxes on hundreds of unique elements 
comprising T/P [Rosborough and Antreasian, 1991]. Once the fluxes on each surface are 
determined for all locations around a single orbit, their transient temperatures were 
determined. These fluxes and temperatures are then converted to resultant forces based 
on the surface properties and orientation of each element. The resultant forces were then 
summed to get the total force acting on the spacecraft at each point in the orbit. A wide 
array of orbit orientations were evaluated for each force and acceleration histories in the 
radial, along-track, and cross-track directions were generated. The direct solar radiation, 
Earth albedo and IR micro-models were evaluated at B’ increments of 4 degrees and 30 
points about the orbit for each B’. The thermal imbalance acceleration on the spacecraft 
was computed at B’ increments of 8 degrees and 72 orbit angles for each B\ 
Consideration also has to be given to the reflective and absorptive properties of the 
surfaces which will degrade with time due to the inhospitable orbit environment. 

A few assumptions have been made in computing the T/P micro-models. These include: 

• Circular orbit 

• Ramp and flip yaw steering modes ignored 

• Constant ‘beginning of mission’ surface properties 

• Diffuse reflection and emission obey Lambert’s cosine law 
for angular distribution of energy 

• The contribution of the force due to the specular 
surface-to-surface reflections is negligible 

• Cylindrical Earth shadow model 

• Constant solar, Earth albedo, and Earth IR fluxes 

• All components of the spacecraft reflect the same percentage of 
reflected energy specularly 

The last assumption is due to data not being currently available for each surface coating, 
but will be incorporated later. Philosophically, the circular orbit and simplistic yaw 
steering provide a useful intermediary "orbit" on which to derive the response of the 
spacecraft to the indicated forces. 

Macro-Models 

The micro-models yield the best available detailed force model history of the spacecraft. 
However, their computational requirements are excessive, taking as much as 1 day to 
process a single orbital revolution on a typical computer workstation. Therefore, a more 
simplistic approach which still mimics the spacecraft’s reaction to radiative forces at a 
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satisfactory level has been adopted. This concept is based on approximating the satellite 
shape with a combination of flat plates. The nonconservative forces acting on each of 
the composite surfaces are computed independently. All plate interaction effects, such 
as shadowing, reflection, and conduction are ignored. This yields vector accelerations 
which are summed to compute the total effect on the spacecraft center-of-mass. The 
algorithm includes the ability to adjust aggregate parameters associated with each flat 
plate to obtain a better fit to the actual satellite acceleration history based on orbit errors 
sensed from laser tracking data and from telemetered satellite on-orbit temperatures. For 
T/P, a box-wing shape was chosen with the plates aligned along the satellite body-fixed 
coordinate system (Figure 1). 

During development, every attempt has been made to keep the macro-models as simple 
as possible while still meeting T/P mission requirements and avoiding telemetry 
dependence. These concepts are often in conflict and, as in the case of the spacecraft 
radiation model, simplicity yielded to the needs of meeting modeling requirements. 


Solar and Earth Radiation 

Solar, albedo, and infrared emissions are the three external radiative fluxes acting on a 
spacecraft. The radiation pressure acting on a flat plate can be computed using the 
following equation [Milani et a!., 1987], assuming a Lambertian diffusion, 

r= _GAcos0 £ 2 ($/3+p COS 0) jj+(i-p) s] 

Me 

where 

r as radiation pressure acceleration on the flat plate 
A a surface area of the flat plate 

e a angle between surface normal and Sun source vectors 

G s radiation flux from source 

n a surface normal vector 

s a source incidence vector 

p a specular reflectivity (percentage of total incoming radiation) 

& a diffusive reflectivity (percentage of total incoming radiation) 

M a satellite mass 

c a speed of light 

The adjustable parameters are area and specular and diffusive reflectivity. These 
parameters are averaged values which represent the consolidation of the spacecraft’s 
complex shapes and material properties into a single, homogeneous flat plate. 
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The albedo and infrared accelerations use a similar acceleration equation as the solar 
radiation. However, the flux magnitude is different. Also, the source vector is the Earth 
grid spot-to-satellite vector rather than the solar incidence vector. The spot definition and 
location are defined by Knocke et al., [ 1988]. Note that this model is not self-shadowing. 
The total albedo/infrared acceleration can be expressed as: 

(2 (vs^cose^vu-p,) 3,1 


where i ■ plate of interest 

j a Earth spot of interest 

n a total number of plates 

Spacecraft Radiation 

Two separate types of fluxes affect the flat plate temperatures: internal and external. 
Internally, the equipment dissipates radiation which serves to heat the satellite surfaces. 
Externally, the solar radiation, albedo, and infrared fluxes cause surface heating. The 
force exerted on a surface due to thermal emission, assuming a Lambertian diffusion 
function, can be expressed as: 


F=- 


2 Ao 
3c 


eT 4 fi 


where 


c 3 emissivity 

a a Stefan-Boltzmann constant (5.67E-08 W/m 2 /K 4 ) 
A a surface area 

T a temperature (Kelvin) 

c a speed of light (m/s) 

n a surface normal vector 


The temperature history algorithm, however, is not as clearly defined. One must take into 
account the complexities of (1 )occultation effects, (2)oblique illumination, and (3) the 
spacecraft’s thermal inertia, without losing sight of the need for simplicity and generality. 


The temperature for a surface exposed to sunlight is modelled as 

! (l) H't). 


T = a + c*cos 
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and while in shadow as 


T = a + 



t 2 +s 2 


where 


s a - 




a a cold equilibrium surface temperature 

c a delta temp, between cold and hot equilibrium 

d a transition time from hot to cold equilibrium temp. 

f s transition time from cold to hot equilibrium temp. 

x a rotation rate/thermal inertia constant 

f, a time since shadow exit 

t 2 a time since shadow entry 

s f a shift parameter to ensure continuity 

s 2 a shift parameter to ensure continuity 

9 — angle between surface normal and solar incidence 

s angle between normal and Sun vectors at shadow entry 

9 ^ = angle between normal and Sun vectors at shadow exit 


The adjustable parameters are area, emissivity, and all five temperature terms ( a,c,d,fji ). 

An explanation of these equations is appropriate. First, note that solar radiation is the 
only direct effect influencing the temperature. That is to say, the 9 angle and the time 
parameters are based only on solar illumination and neglect the albedo and IR effects. 
However, albedo and IR do indirectly influence the apriori values of all the adjustable 
temperature parameters. A plate’s orientation with respect to the Sun dictates which 
temperature algorithm to use. The cosine term in the sunlight equation allows for the fact 
that an obliquely illuminated plate will have a lower temperature than one perpendicular 
to the sun vector. This is illustrated in Figure 3 which shows the temperature profile at 
B’= 0° for a node on the X+ plate as predicted by the micro-model. At an orbit angle 
of 330°, the face enters sunlight and warms rapidly, reaching maximum exposure at Q = 
0°. The cooling from this point until the plate enters shadow at fi = 90° results from the 
decreasing projected area as it rotates with respect to the sun. In contrast, a plate’s 
cooling pattern, when occulted by either the Earth or the spacecraft, is independent of 
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the sun position and, therefore, no such allowance needs to be made. The x parameter 
in the denominator of the cosine term accounts for the fact that the temperature is not 
directly proportional to the rate at which spacecraft rotation moves a plate from direct to 
oblique illumination. Without x, there is no delay in reaching the cold equilibrium value as 
the plate enters shadow. For example, Figure 3 shows that at shadow entry, n = 90°, the 
plate’s temperature is not at its coldest point. 

The exponential term addresses the occultation transition effects. As a face enters or 
leaves shadow, its temperature can be approximated by an exponential curve as shown 
in Figure 4, depicting the B’= 0° micro-model temperature profile for a node on the X- 
face. A different time constant (d or f) is applied depending on whether the surface is 
heating or cooling. 

Finally, a shift term is introduced to ensure continuity in the transition from the sunlight 
to shadow temperature equation. This assumes that the plate will reach its cold 
equilibrium temperature in shadow before heating begins. Given that this assumption is 
not true, a different set of shift parameters must be established. For example, during 
sinusoidal yaw the X- face of T/P is occulted only by the Earth. As B’ increases, this 
shadow time gradually decreases. Therefore, the plate will always reach its hot 
equilibrium temperature and not necessarily reach its cold equilibrium temperature as 
exhibited in Figure 5. In this case the following temperature history algorithms are used: 
Sunlight: 


T = a + 


c * cos/— 
\x 




where 


Shadow: 


where 
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More simplistic representations have been tried with varying success. However, they all 
fail to replicate certain thermal behaviors exhibited by the micro-model. A model of this 
complexity is necessary to meet T/P mission requirements and, therefore, this is the 
chosen parameterization for T/P modeling. 


Results 

In order to test the validity of these macro-models, a comparison has been performed 
between the acceleration histories predicted by the macro and micro-models. A Bayesian 
least squares estimation procedure has been used to tailor the adjustable macro-model 
parameters to better fit the micro-model generated acceleration histories for the solar 
radiation and the thermal imbalance nonconservative forces as outlined by Marshall et 
at., [1991]. Apriori values with realistic uncertainties served to constrain some certain 
highly correlated parameters that could not be recovered independently. Specifically, the 
solar array parameters were not adjusted since their properties are relatively well known 
and they do not represent an average of many smaller surfaces of varying characteristics. 
Also, to ensure realistic temperature values, the equilibrium temperatures, a and c, were 
constrained so as not to stray more than a few degrees from the values predicted by the 
micro-model. Parameters associated with the X+, Y+, Y-, and SA- faces exhibited the 
weakest recovery due to their limited solar exposure. To date, the solar radiation and 
thermal imbalance forces have been fit independently and appropriate parameters 
recovered. Therefore, nonphysical properties could result when the terms are considered 
jointly. With the delivery of the "as-built" spacecraft models, new micro-models will be 
generated and a more aggressive and thorough macro-model parameter recovery will be 
undertaken. These adjusted values will be adopted as nominal values in the actual 
precision orbit determination computations. 

All of the macro-models described herein have been implemented in the precision orbit 
determination software package at NASA/GSFC, GEODYN [Putney et at., 1991]. The 
results presented in this paper use GEODYN and the T/P macro-models to simulate a 10 
day T/P orbit. Table 4 describes the individual plate normal vectors in the satellite body- 
fixed system used in the "box-wing" representation of the T/P spacecraft. Tables 5 and 
6 show the actual macro-model parameter values used in this testing as derived from the 
least- squares adjustment outlined previously. It is important to note that a constant value 
of .34 for Earth albedo and .68 for Earth emissivity were used to represent the Earth 
albedo and IR forces. These values were used as constants to be consistent with the 
micro-model generation. 
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TABLE 4. Macro-Model Plate Normal Vectors in the Spacecraft Body-Fixed System 


PLATE 

X 

Y 

z I 

X+ 

1.0 

0.0 

0.0 

X- 

-1.0 

0.0 

0.0 

Y+ 

0.0 

1.0 

0.0 

Y- 

0.0 

-1.0 

0.0 

Z+ 

0.0 

0.0 

1.0 

Z- 

0.0 

0.0 

-1.0 

SA+ 

1.0 

0.0 

0.0 

SA- 

-1.0 

o 

o 

0.0 


Table 5. Plate Optical and Thermal Characteristics for Tuned Thermal Model 



X+ 

mm 

Y+ 

di 

Z+ 

Z- 

SA+ 

> 

1 

| Area 

4.71 

4.71 

8.18 

8.18 

8.32 

8.32 

25.50 

25.50 I 

| Specular Ref. 

.201 

.244 

.886 

.782 

.239 

.275 

0.05 

0.17 

| Diffuse Ref. 

.375 

.386 

.302 

.339 

.390 

.363 

0.22 

0.66 

Emissivity 

.769 

.995 

.873 

.714 

.770 

.746 

0.87 

0.88 

Temp. A 

181 

168 

191 

190 

240 

103 

236 

234 

Temp. C 

233 

178 

18 

63 

98.5 

125 

110 

96 

Time D 

621 

282 

759 

426 

519 

680 

805 

806 

Time F 

111 

120 

624 

487 

767 

413 

828 

866 

ThetaX 

1.25 

1.00 

1.05 

1.00 

1.06 

1.15 

1.00 

1.00 
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Table 6. Plate Optical and Thermal Characteristics for Tuned Solar Radiation Pressure 

Model. 


1 

X+ 

13 

Y+ 

■3 

Z+ 

Z- 

SA+ 

SA- | 

Area 

3.74 

3.77 

8.27 

8.07 

8.67 

8.44 

21.4 

21,44 i 

Specular Ref. 

.201 

.244 

.886 

.782 

.239 

.275 

0.05 

0,17 1 

Diffuse Ref. 

.375 

.386 

.302 

.339 

.390 

.363 

0.22 

0.66 1 

1 Emissivity 

.769 

.995 

.873 

.714 

.770 

.746 

0.87 

0.88 

1 Temp. A 

181 

168 

191 1 

190 

240 

103 

236 

234 

1 Temp. C 

233 

178 

18 

63 

98.5 

125 

110 

96 

Time D 

621 

282 

759 

426 

519 

680 

805 

806 

Time F 

111 

120 

624 

487 

767 

413 

828 

866 

ThetaX 

1.25 

1.00 

1.05 

1.00 

1.06 

1.15 

1.00 

1.00 


The most important of all the T/P specific models added to the orbit determination 
software is the attitude control logic because all the subsequent nonconservative force 
models rely on the correct orientation of the individual plates. Therefore, substantial effort 
was expended checking this logic. The results not only validate the code but provide 
penetrating insight into the T/P force modeling. 

Figures 6-8 demonstrate the behavior of the B\ yaw, and solar array pitch angles for a 
spacecraft in a simulated T/P orbit over a 1 -year period. These parameters were sampled 
at a rate considerably larger than the orbital period. Therefore, these plots do not 
attempt to show physical dependence on orbital period or ramping region data. Instead 
they demonstrate parameter amplitude and trends. Figure 6 exhibits the variations in B’ 
over 1 year. During this period, B’ variations complete a full cycle. The absolute value 
of B’ never exceeds 89.5° since the orbital inclination is 66° and the ecliptic angle is 23.5°. 
Figure 7 demonstrates that all the yaw steering regimes are exercised and that the 
model adheres to the logic specified in Table 3. Note that the yaw angle transitions 
directly from 0° to -180° instead of following the dictated flip algorithm as the spacecraft 
crosses B’=0°. It is assumed that tracking data taken during this regime is ignored and 
that its effect on the nonconservative force modeling would be minimal. In Figure 8, note 
that the solar array pitch angle (SAPA) oscillates about 0° or 1 80° depending upon the 
yaw steering algorithm. In fixed yaw, where the SAPA is directly related to the orbit angle, 
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one sees the largest variations in amplitude. In sinusoidal yaw, where the intent is to 
keep the solar array pointing perpendicular to the Sun, the SAPA amplitude gradually 
decreases with increasing B\ Almost no variation in SAPA is seen at high B’. 

Figures 9-1 5 show the projected area of each face in the Sun direction as a function of 
B’ and orbit angle. These plots, more than anything else, lend to the understanding of the 
complexities of the T/P radiative nonconservative forces. Beginning with the X+ face in 
Figure 9, this face is illuminated as the spacecraft leaves the Earth occupation boundary 
until it reaches a 90° orbit angle. Once the satellite enters sinusoidal yaw, this plate is 
never again illuminated. This is done to keep the propulsion module, or X- face, in the 
sun direction and, therefore, keep the propellants warm, for as much of the time as 
possible, as shown in Figure 10. Note that the Earth occupation effect gradually 
decreases with increasing B\ At B’>56°, the Earth’s shadow no longer occults the T/P 
orbP. Figures 1 1 and 12 demonstrate how little the Y+ and Y- faces are illuminated. The 
Y- face receives more sunlight on average than the Y+ as P is on the same side of the 
spacecraft as the solar array. Figure 13 exhibits the Z+ plate exposure. This nadir- 
pointing face is unique in that for B’<56° P is illuminated twice per orbit. From 0° to 180° 
orbP angle the face is rotated away from the sun while the spacecraft itself is illuminated. 
It then receives sunlight for the brief periods between Earth occupation and these 180° 
and 360° orbit angles. These periods of sunlight increase wPh increasing B’. Figure 1 4 
shows the very symmetric illumination pattern of the Z- face. Note that, as with the Z+ 
face, the sun’s rays become more glancing as B’ increases and the projected area goes 
to nearly zero. Finally, Figure 15, showing the solar array projected area, clearly outlines 
the Earth occupation region. P is the largest and most stable, in terms of projected area, 
of all the plates. These visibility plots are the main drivers behind the radiative force 
model acceleration histories discussed in the following sections. 


Solar Radiation 

Solar radiation is the dominant radiative force acting on T/P. Ps acceleration profile is 
driven by the large area of the solar array which tracks the Sun throughout the orbP. 
Plots of the micro-model solar radiation acceleration in the along-track, cross-track and 
radial directions are shown in Figures 16, 18, and 20. Data spacing is at 4-degree B’ 
increments and 15 degree orbP angle increments. The "U" shaped outline of the Earth 
occupation region is clearly visible. At B' = 0° the spacecraft is in fixed yaw and the Sun 
is edge on to the orbP plane. Therefore, the plots correctly display no cross-track 
acceleration, a sinusoidal signal in the radial direction, and a similar signal offset 90 
degrees in phase in the along-track component. T/P is in a sinusoidal yaw at B' = 4(F 
and the almost constant cross-track force drops to zero during occupation. The along- 
track and radial accelerations have smaller amplitudes but virtually the same shape as at 
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B' = (f. At 8' = 88°, the spacecraft is in continuous sunlight, the sun is nearly 
perpendicular to the orbit plane, and the cross-track acceleration dominates. 

As exhibited in Figures 17, 19, and 21, the macro-model was exercised in GEODYN to 
generate data for comparison with the micro-model. Notice that the macro-model 
captures all of the features of the micro-model, even though the profiles change 
drastically throughout the different B’ regimes. The residuals between the micro and 
macro-models in each of the three directions are displayed in Figures 22-24. The 
prominent residual spikes appearing around the Earth occupation boundary are caused 
by a different definition of this boundary’s location in the micro and macro-models and, 
therefore, is not a problem since GEODYN calculates shadow entry and exit times very 
precisely. It is anticipated that the remaining residual exhibited in the along-track and 
radial residual plots will be greatly reduced with proper tuning of the solar array macro- 
model parameters. 


Farth Albedo 

The size of the Earth albedo acceleration on the T/P spacecraft is an order of magnitude 
smaller than that of the solar radiation pressure. Its signal is quite different from the solar 
radiation in that the solar array is not always oriented nearly perpendicular to the incident 
flux vector. The micro-model albedo acceleration history in the along-track, cross-track, 
and radial directions is shown in Figures 25, 27, and 29. The accompanying macro- 
model and residual plots are shown in Figures 26,28,30 and 31-33, respectively. Note 
that the acceleration is dominated by the radial direction peak at 90° orbit angle and 0° 
B’ (Figures 29,30). At this point the maximum spacecraft projected area (solar array and 
Z+ face) is perpendicular to the albedo flux. This signal gradually fades as the orbit 
angle approaches 0° or 180° or as B’ increases, and the solar array moves edge on to 
the incident flux. Similarly, the along-track acceleration has the greatest amplitude at low 
B*. The opposing humps in Figures 25 and 26 represent similar spacecraft maximum 
projected area configurations on either side of the 90° orbit angle boundary. From 0° to 
90°, the albedo pushes the spacecraft, while from 90° to 180° it slows the spacecraft 
motion (Figures 25,26). The residual plots do not exhibit the same spiking at the 
occultation boundaries as seen in the solar radiation plots since there is virtually no 
albedo flux striking the spacecraft in this regime. It is important to remember two points 
when examining these residual plots. First, the micro-models were generated using a 
B’=10° fixed to sinusoidal yaw transition boundary. The macro-models, however , used 
a B’=15° boundary value, this leads to some of the signal seen in the residuals for this 
region. Second, the albedo and IR macro-models were never tuned to the micro-models. 
Therefore, a much better fit is anticipated in the future. 
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Earth IR 

The Earth IR acceleration is very similar to the albedo signal. However, it is not influenced 
by Earth occultation and has nearly constant flux throughout the orbit. Therefore, the 
accelerations exhibit a symmetry about 1 80° orbit angle. The micro-model, macro-model, 
and residual acceleration plots are displayed in Figures 34-42. The narrow lip or shelf at 
low B’ in the cross-track plots exists because, until the spacecraft transitions to sinusoidal 
yaw, there is virtually no cross-track signal. Also, observe that the largest magnitude 
along-track accelerations occur at low orbit angles in the macro-model (Figure 35) but 
at high orbit angles in the micro-model (Figure 34). The dips in the radial acceleration 
in the micro-model (Figure 38) are symmetric whereas they are not symmetric in the 
corresponding macro-model plot (Figure 39). This is because the macro-model assumes 
that the reflective characteristics of the material behave in the same way in both the visible 
and infrared wavelengths. This assumption is not correct and will be addressed in the 
future study of the "as-built" micro-models. Finally, be aware that the same problems with 
the different yaw steering boundaries and the absence of tunes parameters discussed 
for the albedo also apply for the IR force. 


Spacecraft Radiation 

The thermal imbalance acceleration on the T/P spacecraft is the most difficult of the 
nonconservative forces to model. However, using the algorithms defined previously, the 
macro-model approximates the micro-model acceleration profiles very well as shown in 
Figures 43-51. The micro-model plots were generated with data spaced at 8 degrees 
in B’ and 5 degrees in orbit angle. The overall signature of this force is similar to the 
solar radiation pressure, which is expected since both forces are primarily driven by the 
spacecraft-Sun orientation. However, Earth albedo and IR have some effect, especially 
at low B’ where they have the greatest magnitude, as demonstrated in temperature plots 
(Figures 52-67). 

The large residuals between the micro and macro-models in the along-track and radial 
directions (Figures 49 and 51) at low B’ are attributed to the fact that the thermal 
imbalance macro-model does not consider heating due to Earth albedo or IR. This signal 
dies off quickly as the spacecraft moves into sinusoidal yaw where the solar array has 
less projected area in the planet direction and the X- face is oriented toward the sun at 
all times. Also, the discontinuities seen in the along-track macro-model (Figure 44) at this 
yaw steering boundary are due to instantaneous changes in the macro-model plate 
shadow boundary location. This is most pronounced for the X- face in the 0° to 90° orbit 
angle region where the plate goes from shadow to maximum illumination as soon as the 
spacecraft crosses B’= 15°. This same discontinuity effect, coupled with differences in the 
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shadow boundary location discussed previously, acts on the Z+ plate to create the large 
spike near B'=56° in the radial residual plot. 

As a further check of the model, the predicted temperature profiles for each plate were 
compared to selected nodes from the micro-model. Recall that the spacecraft radiation 
acceleration is driven by the temperature gradient across the plates. Reasonable 
accelerations can be obtained with unreasonable absolute temperatures. By looking at 
the temperature behavior emulated by the tailored model, the physical reasonableness 
of the resulting model is assessed. These comparisons are shown in the temperature 
profiles shown in Figures 52-67. Note that the micro-model temperature profiles are from 
select nodes on the spacecraft and may not be representative of what is happening on 
the entire plate. Figures 52-55 demonstrate the macro-model’s excellent ability to predict 
the solar array temperature in all regimes. This is attributed to the fact that the solar 
array’s material properties are nearly homogeneous, it has a low thermal inertia, and it 
is oriented for maximum solar exposure. The X+ face is shown in Figure 56 and 57. 
Notice that it experiences wild fluctuations at B’=(f but quickly reaches an almost 
constant cold equilibrium value as the spacecraft enters sinusoidal yaw and the plate is 
no longer illuminated. The warming due to Earth IR can be seen between approximately 
235° and 325° n at low B’ on the micro-model plot only (Figure 56). This is also seen in 
the Z- plots (Figures 66-67). As discussed previously .this effect is not duplicated in the 
macro-model. Figures 58 and 59 show the X- temperature profile which exhibits the 
same large temperature fluctuations at low B’ as experienced by the X+ face. Note the 
significant change in the profile as soon as the spacecraft begins sinusoidal yaw steering 
(j 3 ‘= 15 ?y This has a pronounced effect on the thermal imbalance acceleration as 
outlined previously. Yaw steering causes the dip centered at n=90° while Earth 
occultation causes the larger temperature drop between n=23S° and n=325 ° . Both of 
these effects disappear as B’ increases above 56° where the satellite experiences 
continuous sunlight and has very little yaw motion. The Y+ temperature profile shows 
significant deviation between the micro and macro-models. The macro-model values 
make physical sense when the plate’s solar illumination is considered. The temperature 
consistency shown in the micro-model suggests this particular node is not representative 
of the Y+ plate and that it might be located in or near the spacecraft interior, where the 
temperature is known to be nearly constant. The Y- plate also demonstrates difference 
between the two models, but to a lesser degree. This face, located on the same side as 
the solar array, experiences illumination and, therefore, heating from rays reflected from 
the solar array. This interaction is neglected in the macro-model. The Z+ face (Figures 
64-65) best demonstrate the albedo effect on the temperature signal in the n=0° to 180 f 
range. The other half of the orbit shows the rapid heating and cooling associated with 
the plate illumination twice per orbit. Overall, the macro-model performs more than 
adequately in both acceleration and temperature space. Further verification will be 
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possible on-orbit where direct comparisons with telemetered temperature data can be 
performed. 


Modelabilitv and Error Analysis 

The results presented in the previous sections have given a measure of the macro-model 
success in modeling satellite accelerations. Mission requirements, however, require that 
model performance be evaluated in terms of radial orbit error. In order to quantify the 
radial orbit error produced by macro model errors, the following analysis was performed. 

The micro-model acceleration histories are considered by default to be "truth". Certainly, 
they are the result of a very rigorous analysis and are the best representation of "truth" 
that is currently available. Thus, any mismodeling by the macro-model is represented by 
the macro-micro model acceleration residual histories previously discussed. During this 
analysis, each of the macro model force errors were considered individually. This was 
necessary, since each force was individually tuned to the micro-model. The exception 
is the Earth albedo and IR, which simply use the solar model tuned parameter values 
found in Table 5. 

GEODYN was modified to include a routine that does a bilinear interpolation (over B' and 
ft) on the micro-model acceleration data for each of the forces considered (solar, albedo, 
IR, thermal). Thus, micro-model accelerations were computed at integration step time 
within the GEODYN software. Figures 68 to 79 show the bilinearly interpolated micro 
model acceleration histories for a B' region of 10.5 to 36 degrees. Comparisons of these 
figures with the micro-model acceleration histories discussed in previous sections show 
that the bilinear interpolation is adequate. The "micro-model GEODYN" version produced 
four data sets containing true-of-date X, Y, Z orbit data for a 10 day T/P arc, using each 
of the micro-model accelerations in turn. A 30 second integration step size was used for 
all runs, and the data generation interval was 180 seconds. The particular 10 day arc 
spans a B' region of 10.5 to 36 degrees. Therefore, this arc covers the fixed yaw to 
sinusoidal yaw regime, and the analysis is weighted towards the acceleration residuals 
in this region. Four separate data reduction runs using the macro-model accelerations 
were made on the four micro-model generated data sets. Thus, for each individual force 
considered, the only difference in the force modeling between the data generation and 
data reduction runs is the particular macro-micro model difference. Only the state, a 
single drag coefficient, and a single solar radiation pressure coefficient were adjusted 
over the 1 0-day arc interval. The adjustment of specific macro-model parameters were 
not part of this study, and will be addressed extensively in a companion report in 
preparation. The radial RMS orbit error was computed from the residuals of the macro- 
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model-fit orbit to the micro-model-generated orbit. 

Table 7 gives the modelability and error analysis summary. TTie RMS radial orbit errors 
over a 10 day T/P arc for each of the individual macro-model forces is less then 6 cm, 
demonstrating that each individual macro-model meets mission requirements. It should 
be stressed that these results were achieved without adjusting T/P macro-model specific 
parameters, and that the inclusion of these parameters should improve modelability and 
reduce the orbit errors. 

Some comments about the results in Table 7 are given below. The modelability of the 
along-track and cross-track albedo and IR accelerations is poor. However, the surface 
properties used to describe the albedo and IR models were never tuned to the micro- 
models, and currently the solar model parameter set is employed. The modelability of 
the accelerations should improve when macro-model parameters are adjusted. 
Furthermore, the radial component of the albedo and IR accelerations is an order of 
magnitude larger than the other components. The modelability of this component is quite 
good. The along-track thermal acceleration also demonstrates deficiencies due to 
problems in modeling and tuning the X- plate temperature. Further model tuning is 
necessary for this X- face. The solar radiation pressure is currently modeled at the 10% 
level and is anticipated to improve with the tuning of the solar array diffuse and specular 
reflectivity parameters. 
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Table 7. Modelabllity and Error Analysis Summary 


FORCE 

RMS of Micro 
Model Force for 
B' of 0° to 88° 
and Orbit Angle 
of 0° to 360° 

RMS of Macro- 
Micro Model 
Residuals 

RMS Radial 10 
Day Orbit Error 
(For All 
Components) 

Solar Along-Track 
Solar Cross-Track 
| Solar Radial 

3.2x10* M/S 2 
4.5x10* M/S 2 
2.7x10* M/S 2 

4.3x10* M/S 2 
4.3x10* M/S 2 
4.9x10* M/S 2 


I Total Radial Orbit Error for 1 0 Day Arc -* 


5.1 cm 

1 Albedo Along-Track 
| Albedo Cross-Track 
| Albedo Radial 

6.1x10 10 M/S 2 
8.3x10 10 M/S 2 
4.6x10* M/S 2 

2.3x1 0 10 M/S 2 
3.3x1 0 10 M/S 2 
2.4x1 O' 10 M/S 2 


| Total Radial Orbit Error for 10 Day Arc -* 


2.1 cm 

| IR Along-Track 
| IR Cross-Track 
| IR Radial 

6. 1x1 0 10 M/S 2 
7.4x1 0 10 M/S 2 
5.6x10* M/S 2 

5.4x1 0’° M/S 2 
5.3x10 10 M/S 2 
5.4x10 ’° M/S 2 


1 Total Radial Orbit Error for 1 0 Day Arc -*• 


2.2 cm 

Thermal Along-Track 
Thermal Cross-Track 
Thermal Radial 

2.0x10* M/S 2 
3.6x10* M/S 2 
1.6x10* M/S 2 

1.0x10* M/S 2 
5.9x1 0 ’° M/S 2 
5.5x1 0 ’° M/S 2 


| Total Radial Orbit Error for 10 Day Arc -* 


5.2 cm 


The above analysis shows that with no adjustment of the macro-model surface-specific 
parameters and a minimal adjusted parameter set, the individual macro-model errors 
meet mission requirements. However, the error produced by the sum of the 
nonconservative forces must meet the 6 cm radial RMS mission requirement. This 
analysis could not be properly performed since a combined force parameter set has not 
yet been developed but will be addressed in a companion paper. In order to get an 
estimate of the combined force error, the above analysis was performed using the 
combined micro and macro-model forces. The thermal model tuned parameter set was 
used in this analysis. Adjusting the state, a single drag coefficient, and a single solar 
radiation pressure coefficient the radial RMS orbit error of a 10 day arc was 6.98 cm. A 
short run was made including the adjustment of the solar array front diffuse reflectivity 
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parameter and demonstrated marked improvement, giving a flavor for the type of 
improvement expected when macro-model parameters are adjusted. 
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CONCLUSIONS 


Summary 

In order to meet precision orbit determination requirements for geodetic satellite missions, 
specifically TOPEX/Poseidon, detailed models of the radiative forces acting on the 
spacecraft have been constructed. Solar, Earth albedo. Earth infrared, and the 
spacecraft’s thermal radiation effects have all been considered. A detailed finite element 
analysis has been performed to compute the total force and induced accelerations acting 
on the satellite. This required a precise description of the satellite shape, material 
properties, and attitude control algorithm. Because these models are too computationally 
intensive to be incorporated into the orbit determination software, a more simplistic model 
which approximates the finite element acceleration profiles has been developed. It is 
based on depicting the satellite as a combination of flat plates and computing the 
nonconservative forces acting on each plate independently. These acceleration vectors 
are summed to produce the overall effect on the satellite center-of-mass. For T/P, a box- 
wing shape is used. Each plate has associated parameters which can be adjusted to 
improve model performance with respect to the micro-model analysis and, when the 
spacecraft is on orbit, to the tracking observations. The adequacy of these macro- 
models have been assessed through direct comparison with the micro-models. The 
effect of these residuals on the orbit error budget has been studied. Analyses indicate 
these precise models individually meet mission requirements. It has been shown that 
improved modelability will be achieved when macro model specific parameters are 
adjusted. 


Future Work 


With the incorporation of the macro-models into the GEODYN orbit determination code 
[Putney, 1991 J, an intensive study of parameter recoverability and correlation can and will 
be undertaken. Once Fairchild supplies the T/P "as-built" models, all of the micro-models 
will be reiterated and the box-wing parameters will be tailored to these reiterated 
solutions. After launch the actual on-orbit satellite performance as sensed from l as er 
tracking and telemetered temperatures will be used to further refine these values. 

Additional error analysis with the "as-built" tuned macro-model parameters needs to be 
performed. This will enable the determination of the proper orbit error due to the sum 
of the macro-model force errors. Macro-model specific parameters need to be adjusted 
to improve modelability. Finally, error analysis using simulated laser tracking must be 
performed. This is because the model performance is to be evaluated in terms of radial 
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orbit error in computations supported exclusively by laser tracking. 

When adding the macro-models to the GEODYN code, every attempt is being made to 
make the code as general as possible rather than T/P specific. Any satellite specific code 
is isolated to facilitate the change to other spacecraft. Therefore, these same algorithms 
can be used to support improved orbit determination for other satellites. Specifically, two 
European satellites, SPOT-2 and ERS-1 , are well suited to this task as are the satellites 
forming the United States’ GPS constellation. 

With improvements in static gravitational and tidal force models, our attention must now 
focus on other contributing orbit error sources. The days of treating complex geometric 
satellites as cannonballs in orbit determination are over. The model development 
described in this paper represents a first step in addressing what may become the most 
dominant error source in precision orbit determination in the near future. 
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Figure 1 . (a) The TOPEX/Poseidon Spacecraft, (b) Micro-Model Approximation, 

(c) Macro-Model Approximation 
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Figure 2. TOPEX/Poseidon Inertial Coordinate System (Perrygo, 1987) 


MICRO-MODEL TEMPERATURE PROFILE 



Figure 3. T/P X+ Plate Node Temperature Profile Over One Revolution at B'=0° 
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MICRO-MODEL TEMPERATURE PROFILE 



Figure 4. T/O X- Plate Node Temperature Profile Over One Revolution at B'=0° 


MICRO-MODEL TEMPERATURE PROFILE 
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Figure 5. T/P X- Plate Node Temperature Profile Over When Revolution at B' = 40° 
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BETA PRIME 



Figure 6<a). B' Variation Over the First 3 Months of a 1 Year T/P Orbtt 
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Figure 6(b). B' Variation Over the Second 3 Months of a 1 Year T/P Orbit 
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BETA PRIME 


TWBO 3 MONTHS 



Figure 6(c). B Variation Over the Third 3 Months of a 1 Year T/P Orbit 


BETA PRIME 



Figure 6(d). B' Variation Over the Fourth 3 Months of a I Year T/P Orbit 
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Figure 7(a). First 3 Months of T/P Yaw Angle Variation Over One Year 
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Figure 7(b). Second 3 Months of T/P Yaw Angle Variation Over One Year 
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Figure 8(a). First 3 Months of T/P Solar Array Pitch Angle Variation Over 1 Year 


SOLAR ARRAY PITCH ANGLE 


SEOOW5 3 MONTHS 


(Thousands) 

MJO 


Figure 8(b). Second 3 Months of T/P Solar Array Pitch Angie Variation 

Over 1 Year 
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Figure fl(c). Third 3 Months of T/P Solar Array Pitch Angie Variation Over 1 Year 
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Figure 8(d). Fourth 3 Months of T/P Solar Array Pitch Angle Variation Over 1 Year 
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Figure 16. Micro-Model Solar Radiation Pressure Along Track Acceleration 



Figure 17. Macro-Model Solar Radiation Pressure Along Track Acceleration 
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Figure 18. Micro-Model Solar Radiation Pressure Cross Track Acceleration 



Figure 19. Macro-Model Solar Radiation Pressure Cross Track Acceleration 
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Figure 22. Along Track Solar Radiation Pressure Acceleration Residuals. 
(Macro-Micro) Mean = -3.8E-1 1; RMS » 4.3E-9 



Figure 23. Cross Track Solar Radiation Pressure Acceleration Residuals. 
(Macro-Micro) Mean = 1.7E-9; RMS = 4.3E-9 
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Figure 24. Radial Solar Radiation Pressure Acceleration Residuals. 
(Macro-Micro) Mean = 1 .3E-9; RMS = 4.9E-9 
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Albedo Cross Track Acceleration 
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Figure 31. Along Track Earth Albedo Acceleration Residuals. (Macro-Micro) 



Figure 32. Cross Track Earth Albedo Acceleration Residuals. (Macro-Micro) 

Mean = 1.9E-10; RMS=3.3E-10 
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Figure 33. Radial Earth Albedo Acceleration Residuals. (Macro-Micro) 

Mean = 2.2E-10; RMS=2.4E-10 
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Figure 39. Macro-Model Earth IR Radial Acceleration 
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Figure 40. Along Track Earth IR Acceleration Residuals.(Macro-Micro) 

Mean = 8.5E-12; RMS = 5.4E-10 
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Figure 41. Cross Track Earth IR Acceleration Residuals. (Macro-Micro) 

Mean = 4.4E-10; RMS = 5.3E-10 
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Figure 42. Radial Earth IR Acceleration Residuals. (Macro-Micro) 
Mean =- 2.8E-10; RMS = 5.4E-10 
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Figure 43. Micro-Model Thermal Imbalance Along Track Acceleration 



Figure 44. Macro-Model Thermal Imbalance Along Track Acceleration 
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Figure 45. Micro-Model Thermal Imbalance Cross Track Acceleration 



Figure 46. Macro-Model Thermal Imbalance Cross Track Acceleration 
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Figure 47. Micro-Model Thermal Imbalance Radial Acceleration 





Figure 48. Macro-Model Thermal Imbalance Radial Acceleration 
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Figure 49, Along Track Thermal Imbalance Acceleration Residuals. 
(Macro-Micro) Mean = 4.9E-11; RMS = 1.09E-9 
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Figure 50. Cross Track Thermal Imbalance Acceleration Residuals. 
(Macro-Micro) Mean = 5.9E-10; RMS = 7.6E-10 
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Figure 51. Radial Thermal Imbalance Acceleration Residuals. 
(Macro-Micro) Mean = 1.6E-10; RMS = 5.5E-10 
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Figure 64. Z+ Micro-Model Temperature Profile 



Figure 65. Z+ Macro-Model Temperature Profile 





Figure 68. Bilinear interpolated Micro Model Solar Radiation Pressure Along. 

Track Acceleration. For 10.5° to 36° B' 
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Figure 69. Bilinear Interpolated Micro Model Solar Radiation Pressure Cross. 

Track Acceleration. For 10.5° to 36° B' 
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Figure 70. Bilinear Interpolated Micro Model Solar Radiation Pressure 
Radial Acceleration For 10.5° to 36° B' 



Fiqure 71 Bilinear Interpolated Micro Model Albedo Along Track Acceleration. 

For 10.5° to 36° B' 
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Figure 72. Bilinear interpolated Micro Model Albedo Cross Track Acceleration. 

For 10.5° to 36° B' 



Figure 73. Bilinear Interpolated Micro Model Albedo Radial Acceleration. 

For 10.5° to 36° B' 


72 




Figure 74. Bilinear Interpolated Micro Model IR Along Track Acceleration. 

For 10.5° to 36° B' 



Figure 75. Bilinear Interpolated Micro Model IR Cross Track Acceleration. 

For 10.5° to 36° B' 
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Figure 76. Bilinear Interpolated Micro Model IR Radial Acceleration. 

For 10.5° to 36° B' 



Figure 77. Bilinear Interpolated Micro Model Thermal Along Track Acceleration. 

For 10.5° to 36° 6' 
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Figure 78. Bilinear interpolated Micro Model Thermal Cross Track Acceleration. 

For 10.5° to 36° B' 



Figure 79. Bilinear Interpolated Micro Model Thermal Radial Acceleration. 

For 10.5° to 36° B' 
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